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DYNAMICAL INVESTIGATION OF MINOR RESONANCES FOR 

SPACE DEBRIS 

ALESSANDRA CELLETTI AND CAT ALIN GALE§ 


Abstract. We study the dynamics of the space debris in regions corresponding to 
minor resonances; precisely, we consider the resonances 3:1, 3:2, 4:1, 4:3, 5:1, 5:2, 5:3, 
5:4, where a j : ^ resonance (with j, £ G Z) means that the periods of revolution of the 
debris and of rotation of the Earth are in the ratio jji. We consider a Hamiltonian 
function describing the effect of the geopotential and we use suitable finite expansions 
of the Hamiltonian for the description of the different resonances. In particular, we 
determine the leading terms which dominate in a specific orbital region, thus limiting our 
computation to very few harmonics. Taking advantage from the pendulum-like structure 
associated to each term of the expansion, we are able to determine the amplitude of 
the islands corresponding to the different harmonics. By means of simple mathematical 
formulae, we can predict the occurrence of splitting or overlapping of the resonant 
islands for different values of the parameters. We also find several cases which exhibit 
a transcritical bifurcation as the inclination is varied. 

These results, which are based on a careful mathematical analysis of the Hamiltonian 
expansion, are confirmed by a numerical study of the dynamical behavior obtained by 
computing the so-called Fast Laypunov Indicators. 

Since the Hamiltonian approach includes just the effect of the geopotential, we vali¬ 
date our results by performing a numerical integration in Cartesian variables of a more 
complete model including the gravitational attraction of Sun and Moon, as well as the 
solar radiation pressure. 


1. Introduction 

The dynamics of objects moving in the space surrounding the Earth is a subject of 
strong interest, due to the many satellites that have been placed in orbit around our 
planet and that generated many debris (see, e.g., [15]). In this work we are interested to 
the dynamics corresponding to a resonant motion, which occurs whenever the period of 
revolution of the celestial object and the period of rotation of the Earth are commensu¬ 
rable. Resonant motions have been widely used to design the orbit of artificial satellites. 
In particular, two resonances have been fully exploited by the geosynchronous and the 
GPS satellites (see, e.g., [1], [8], [16], [18], [19], [20], [21], [22]). In the first case, the ob¬ 
ject moves in a 1;1 resonance (at about 42164 km from Earth’s center), always pointing 
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the same location on the Earth, since its orbital period is exactly equal to one sidereal 
day. The GPS satellites move in a 2:1 resonance at about 26560 km from Earth’s center, 
which means that they make two orbits during one rotation of the Earth. 

More in general, there is a standard classihcation of different regions of the sky, in terms 
of the altitude above the Earth: LEO (acronym of Low-Earth-Orbit), MEO (Medium- 
Earth-Orbit) and GEO (Geosynchronous-Earth-Orbit) are three regions which divide 
the sky, starting from the Earth’s surface to the geosynchronous ring. Precisely, LEO 
corresponds to the sky between 90 and 2000 km. MEO is the region between 2000 and 
30000 km, which includes GPS as well as other resonances. GEO is located at an altitude 
above 30000 km from the surface of the Earth. 

All three regions are affected by several forces (see Section 2); Erst, the Earth’s at¬ 
traction is very strong and the geopotential must be included within a high degree of 
precision; next, the effect of Sun and Moon is extremely important and must be taken 
into account (see [11]); also the solar radiation pressure plays a relevant role; finally, in 
the LEO region the atmospheric drag must be certainly considered. 

In this paper we concentrate on resonances of lower order (w.r.t. GEO and GPS), 
to which we will refer as minor resonances: 3:1, 3:2, 4:1, 4:3, 5:1, 5:2, 5:3, 5:4, which 
populate the region of the sky between 14000 km and 37000 km from Earth’s center. 
Our study aims at exploiting the minor resonances using mathematical tools based on 
a Hamiltonian approach, which allows us to have a deep understanding of the dynamics 
of such resonances. This task can be accomplished, once we have a model that describes 
with good accuracy the dynamics. To this end, we expand the geopotential to different 
orders, according to the resonance we are considering (see Section 2). However, since the 
expansion might contain a huge number of terms, following [3] we introduce the notion 
of dominant term in a specific region of the orbital elements’ space (see Section 2). This 
allows us to considerably reduce the number of harmonics which really shape the dynam¬ 
ics. For reasonable parameter values, the resonances have a typical pendulum structure, 
showing an island shape surrounding the elliptic point. We present a simple mathemat¬ 
ical algorithm that allows one to compute the amplitudes of the resonant islands with 
a minimum computational effort (compare with Section 3). Gasting together such in¬ 
formation about the dominant terms and the amplitudes of the islands, we are able to 
proceed further in predicting whether the islands associated to the different harmonic 
terms are well separated or they rather overlap giving birth to chaotic motions (the so- 
called splitting or superposition phenomena described in Section 4). The prediction of 
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such behavior is obviously very important, since it could allow for regular or chaotic mo¬ 
tions. Indeed, we also propose a transfer mechanism at low cost, taking advantage of the 
stable or chaotic character of the dynamics as some orbital elements are suitably varied. 
Finally, we study the mechanism of transcritical bifurcations (see Section 5), which occur 
for some resonances and which provoke a sudden change in the stable/unstable behavior 
of the equilibria. 

To summarize, this paper is organized as follows. In Section 2 we present the model, 
using both the Cartesian and Hamiltonian formulations. A measure of the amplitudes 
of the resonances is provided in Section 3. A mechanism of splitting or superposition 
of resonances is given in Section 4, while the occurrence of transcritical bifurcations is 
investigated in Section 5. A model including all main forces, and not just the geopotential, 
is studied in Section 6 using a Cartesian approach. 

2. The model in Cartesian and Hamiltonian formalism 

In this section we introduce the equations of motion of a small body, say S, that we 
identify with a space debris; we assume that S is subject to the influence of the Earth 
and, beside the gravitational interaction, we take into account also the geopotential up 
to a hnite degree. Within the Cartesian formalism, we consider also the effects of Sun 
and Moon, as well as the solar radiation pressure. 


Let us introduce a quasi-inertial frame centered in the Earth. The equation of motion 
in Cartesian coordinates will consider the Earth’s gravitational influence, the geopoten¬ 
tial, the solar attraction, the lunar attraction and the solar radiation pressure. Precisely, 
let us denote by ruE, ms and tum the masses of Earth, Sun and Moon, let Q be the 
gravitational constant. With reference to [3], the equation of motion is given by 


= R3{-e)VV{T)-gms 

• - Pm , rM 


rs 


T - ^5 


+ 


rs 


Gmu 


Ym 


Tm 


+ CrPral (-) 

m |r — r 5 |‘- 


( 2 . 1 ) 


where r, rs, tm represent the position vectors of the debris. Sun and Moon with respect 
to the center of the Earth (see [17] for explicit formulae concerning rs, ^m), R 3 denotes 
the rotation about the polar axis, 6 is the sidereal time, V is the gradient computed with 
respect to the synodic frame, while V{r) is the force function due to the attraction of 
the Earth (see, e.g., [3] for full details). 
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As we can see from ( 2 . 1 ), the contribution of the solar radiation pressure involves the 
reflectivity coefhcient Cr of the debris, the radiation pressure Pr for an object located at 
distance 05 = 1 AU, and the area-to-mass ratio A/m with A the cross-section of the 
debris and m its mass. 

Next we consider just the effect of the geopotential and we provide the corresponding 
Hamiltonian function in terms of the action-angle Delaunay variables (L, G, H, M, cn, D). 
Such coordinates are linked to the orbital elements (a, e, f, M, cu, D), where a is the semi¬ 
major axis, e the eccentricity, i the inclination, M the mean anomaly, uj the argument of 
perigee, D the longitude of the ascending node. Precisely, denoting hj fiE = G^e, one 
has the following relations: 

L = y/ fiEtt , G = Ly/l — , H = Gcosi . 

The Hamiltonian describing the geopotential contribution in (2.1) can be written as (see 

[3]) 

nm G, H, M, UJ, n, 6) = + RearthiL, G, H, M, UJ, n, 6) , 

where the geopotential is given by (see [14]) 

00 n n 00 

Rearth = — , UJ, Q, 6) , (2.2) 

CL ^ a ' 

n =2 m =0 p =0 q=—oo 

where Re is the equatorial radius of the Earth, the well-known inclination and eccentricity 
functions Fnmp, Gnpq are given in [14] through recursive expressions, while Snmpq depends 
on the spherical harmonic coefficients G^^, Snm (see [14]) and on the angle 

'^nmpq = (u - 2p)u + {u - 2p - q)M + m{Vt - 6) . (2.3) 

Let us also introduce the quantities Jnm and Xnm dehned through the relations 

Gnm Jnm eOs(?7lA^n2) , S^m Jnm sin(?T7.Afim) . 

The coefhcients G^m, Snm and Jnm in units of 10“®, as well as the values of A„m, up to 
degree and order 5, are given in Table 1, derived from the EGM2008 model ([10], see also 

[5], [17]). 


In order to provide a description of the resonant motions, we expand Rearth and, 
averaging over the non-resonant terms, we retain the secular and resonant parts, which 
yield the long term variation of the Delaunay variables, hence of the orbital elements. 
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n 

m 

r 

^nm 

^nm 

Jnm 

^nm 

2 

0 

-1082.6261 

0 

1082.6261 

0 

3 

0 

2.53241 

0 

-2.53241 

0 

3 

3 

0.100583 

0.197222 

0.22139 

80°9928 

4 

0 

1.6199 

0 

-1.619331 

0 

4 

3 

0.059215 

-0.012009 

0.060421 

56°1784 

4 

4 

-0.003983 

0.006525 

0.007644 

-14°6491 

5 

4 

-0.0023 

0.000388 

0.00233198 

-2^39321 

5 

5 

0.00043 

-0.00165 

0.001703 

20°9272 

6 

4 

-0.0003256 

-0.0017845 

0.001814 

19°9146 

6 

5 

-0.00022 

-0.00043 

0.000483703 

12°7055 


Table 1. The coefficients Snm, Jnm (in units of 10 and the 
quantities Xnm] values computed from [10]. 


We shall consider the Earth’s gravitational potential up to terms of degree and order 
n = m = N, where N will be given later as it will depend on the specific resonance we 
consider. Let us write Rearth as 

N n n oo 

p _ nsec I pres ^ _ \ ^ \ ^ \ ^ \ ^ 

dearth -^earth -^earth ~ / ^ ^ ^ ^ ' nmpq ? 

n=2 m=0 p=0 q=—oo 

where R^^arthi dearth represent the secular and resonant parts of the Earth’s potential, 
while the coefficients Tnmpq are given in Appendix A. Next we introduce the following 
definition of gravitational resonance. 

Definition 1 . A j : i gravitational resonance for j, i G Z\{0} occurs when the orbital 
period of the debris and the period of rotation of the Earth are commensurable in the 
ratio j/i. In terms of the orbital elements, one has: 

iM-je = 0, jJeN. (2.4) 

Notice that (2.4) is satisfied in concrete astronomical cases within a certain degree of 
approximation and cannot be obviously satisfied exactly. 

By using Kepler’s third law, it follows that a j : £ resonance corresponds to the semima¬ 
jor axis aj:e = {jl£)~^^^ Ogeo, where a^eo = 42164.1696 km represents the semimajor axis 
of the geosynchronous orbit. Table 2 provides the location of the resonances, that we shall 
investigate in this work, as well as those of the 1:1 and 2:1 resonances for comparison. 

For all resonances we write the same expression for the secular part, due to the fact 
that the geopotential coefficient J 2 = J 20 is much larger than any other zonal coefficient 
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J-.i 

a in km 

j-.i 

a in km 

1:1 

42164.2 

4:3 

34805.8 

2:1 

26561.8 

5:1 

14419.9 

3:1 

20270.4 

5:2 

22890.2 

3:2 

32177.3 

5:3 

29994.7 

4:1 

16732.9 

5:4 

36336 


Table 2. Value of the semimajor axis corresponding to several 
resonances. 


(see Table 1): in the expansion of the secular part the most important role is played by 
a term of order 0 {J2)- On the other hand, the resonant parts of the development of the 
geopotential are obtained adding different terms, say Tk for some /c G Z+; we will need 
to compare the strength of snch terms to rednce onr stndy to a fnnction composed by 
the most signihcative contribnting terms, dehned as follows (see [3]). 

Definition 2. Let BJ'g^arth resonant part of Rearth, corresponding to the resonance 

j : i. Let be the associated stroboscopic mean angle. Given the orbital elements 
{a,e,i), we say that a term Tk for some k E of the expansion of , say Tk = 

gk{a, e, i) cos{k + \k) for some function gk and some constant \k, is dominant, if the 
size of \gk{a,e,i)\ is bigger than the size of any other term of the expansion. 

The analysis of the dominant terms allows us to reduce the discussion to a limited 
nnmber of terms as well as to provide an indication of the optimal degree of the ex¬ 
pansions. More precisely, for a given resonance j : i we approximate the Hamiltonian 
function with 

qjresj:i _ - TDsec I r)resj:£ 

^ ' ^earth ’ 

where Rff^rth expanded np to an optimal degree N, which is determined by implement¬ 
ing the algorithm described in [4]. The optimal degree of expansion of Rffadh = J + lj 
except for the resonance 4:1 whose optimal degree is iV = j + 2. The terms which con- 
tribnte to form are listed in Table 3; explicit expressions for the corresponding 

coefficients are given in Appendix A. 

A plot of the dominant terms according to Dehnition 2 for each of the resonances 
considered in this work is provided in Figure 1. 
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N 

terms 

3:1 

4 

7330-2 , 73310 , 73322 , 7431-1, 74321 

3:2 

4 

TsSO-l, 73311, 7430-2, TlsiO, 74322 

4:1 

6 

7441-1, 74421, 7541-2 , 75420 , 75432 , 7642-1, 76431 

4:3 

5 

7440-1, 74411, 7540-2 , 75410 , 75422 

5:1 

6 

7551-2 , 75520 , 75532 , 7652-1, %531 

5:2 

6 

7551-1, 75521, 7651-2 , 76520, %532 

5:3 

6 

7550-2, Tssio, 75522, Tesi-l, 76521 

5:4 

6 

7550-1, Tssil, 7650-2 , 76510 , 76522 


Table 3. Terms whose sum provides the expression of 
order N. 


3. Measuring the amplitude of resonant islands 

In this section we concentrate on the size of the resonant islands associated to the 
dominant terms. First, we introduce in Section 3.1 an elementary mathematical method 
to estimate the size of the resonant island associated to a specihc term, provided that 
we are in a parameter region corresponding to a regular (and not chaotic) behavior (see 
[3]). Examples are given in Sections 3.2 and 3.3. 


3.1. A pendulum-like estimate of the amplitude. Following [3], we sketch an ele¬ 
mentary method which allows us to estimate the amplitude of the island around a given 
j : ^ resonance (see [3] for full details). This estimate is computed by taking into ac¬ 
count the influence of the secular part and just the largest term of the resonant part. 
In what follows, we obtain the width of the resonant island associated to the dominant 
term as a function of eccentricity and inclination. However, it is important to underline 
that in many regions of the phase space (usually for moderate and large eccentricities), 
some resonant harmonic terms with comparable large enough magnitude could coexist. 
Due to a common phenomenon which takes place for almost all minor resonances, called 
splitting of the resonances and detailed in Section 4, these big harmonic terms yield non¬ 
overlapping resonance islands. Therefore, around a given j : i resonance there could be 
multiple resonant islands, according to the values of inclination and eccentricity. 

In this section, we focus our attention on the resonant island having the largest width. 
The resonant Hamiltonian can then be written as 




hi 

'2L2 


+ RZ-,h{L,G,H,uj) + R, 


(3.1) 
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with 

Ni N2 N3 

^ i; 5 : i; R]i’^\L, G, H)cs{h{iM-je) + k2U + k5n) , 

ki=l k2 = l ^ 3=1 

where Ni, N 2 , are integers, denote the Fourier coefficients, cs could be either 

cosine or sine and k = {ki,k 2 ,ks). 

Normalizing the units such that 0 = 1, then from the resonance relation and Kepler’s 
third law, we obtain that the resonant value of the action L is given by 


Lrp.M 


3 


(3.2) 


We expand (3.1) around Lj-es up to second order and we retain only the largest term in 
the resonant Hamiltonian: 

G, H, m - jO, u, Q) = aA- /3A^ + r] cs{k^^%m - jO) + k^'oj + 

with 

A L 

a = a{U,,,G,H,u:)^^ + RT:rtH,L{Lres,G,H,u:) 


res 

2 


^^res ^ 

r] = r]{Lres, G, H) = ’ (3-3) 

where = {k'^°‘^, k^°'^, k^°‘^) denotes the index at which rj is maximum. One can 
show that the variation of A is given by 


AA = 



/S 


From a = /^ie we obtain 


r2 r2 

Aa = — -^ 


k'E 

SO that, using AL = AA, the amplitude of the j : ^ resonant island is given by 

2 Aa=A(^ + 2 L„,./f). 

kE ^ P V P ^ 

with T], (3 as in (3.3) and L^es as in (3.2). 


(3.4) 


We report in Figure 2 the amplitudes of the minor resonances studied in this paper as 
a function of the eccentricity (between 0 and 0.5) and the inclination (between O'’ and 
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Figure 1. Dominant terms (indexes reported in the plots) as a function 
of eccentricity and inclination for the resonances 3:1, 3:2, 4:1, 4:3, 5:1, 

5:2, 5:3, 5:4. 
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90"); in the plots we fixed a; = 0" and = 0°. The color bar provides the size of the 
amplitnde in kilometers. 

In Sections 3.2-3.3 we consider some minor resonances as bench tests for the deter¬ 
mination of the amplitudes using the expression (3.4) and comparing the results with 
an investigation based on the computation of the Fast Lyapunov Indicators (hereafter, 
FLIs), which are defined as the largest Lyapunov characteristic exponents at a finite 
time. FLIs were introduced in [12] and implemented in [3] in the context of space debris 
to which we refer for more details (see also [2] and [4], [13] for cartographic studies based 
on the FLIs). 

3.2. The 3:1 resonance. For the 3:1 resonance, we have five terms defining RH^fh 
Table 3). The amplitude of each dominant term is computed in Table 4 using (3.4) for 
different eccentricities and inclinations. Despite the simplicity of (3.4), the agreement 
with more accurate computations is evident from a comparison with Figure 3 (top and 
middle rows), representing the FLI values as a function of mean anomaly and semimajor 
axis and the bottom row, where the FLI is plotted both as a function of inclination and 
semimajor axis. 

For small eccentricities and small to moderate inclinations, all terms of RUrfh^ except 
Tssio, are small in magnitude, so that a pendulum-like plot is obtained (see Figure 3, 
top left and middle left). The amplitudes of the islands associated to Tasio reported 
in Table 4 are definitely consistent with those computed from Figure 3, top left and 
middle left panels. However, increasing the eccentricity, other terms grow in magnitude 
showing a pendulum structure, although they do not interact with the main resonance 
even for large eccentricities, provided the inclination is small (compare with Figure 3 top 
right). In this case, the estimate (3.4) still provides a good value for the amplitude of 
the resonant island associated to the dominant terms. 

For higher inclinations and larger eccentricities, the main resonance increases a lot in 
amplitude and it interacts with the other resonances, leading to chaotic motions (Figure 3, 
middle right); in this case, as expected, the estimates given by (3.4) do not properly work. 
We notice that the amplitude of the largest term increases significantly in passing from 
i = 10" to i = 30°. In particular, due to the fact that the amplitudes of the main terms 
for i = 10" are not too large and that the center of the different terms are shifted, there is 
no superposition of the resonances (see Figure 3, top right). On the contrary, for i = 30" 
the amplitudes are sufficiently large to provoke an interplay of the resonances generated 
by the different terms (see Figure 3, middle right). This behavior will be the centerpiece 
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Figure 2. The amplitude of the resonances for different values of the 
eccentricity (within 0 and 0.5 on the x axis) and the inclination (within 
0° and 90° on the y axis) for ca = 0°, hi = 0°; the color bar provides the 
measure of the amplitude in kilometers. In order from top left to bottom 
right: 3:1, 3:2, 4:1, 4:3, 5:1, 5:2, 5:3, 5:4. 
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of the discussion of Section 4, where the splitting and superposition of the resonances 
will be explained in detail. 


Dominant term 

e = 0.005, i = 10° 

e = 0.005, f = 30° 

e = 0.5, f = 10° 

e = 0.5, i = 30° 

7330-2 

0.05 

0.05 

5.25 

4.79 

73310 

4.50 

12.57 

5.51 

15.40 

73322 

0 

0.02 

0.23 

1.97 

7431-1 

0.33 

0.46 

3.35 

4.65 

74321 

0.11 

0.52 

1.14 

5.25 


Table 4. Amplitude in kilometers using (3.4) of the dominant terms 
associated to the 3:1 resonances for e = 0.005, 0.5 and i = 10°, 30°. 


The behavior of the amplitude, as computed from the FLI plots, can be obtained from 
the bottom row of Figure 3, which is computed for a hxed eccentricity and a whole 
interval of inclinations (similarly, we could have shown the plots in the (e, a)-plane for a 
hxed inclination). 


3.3. Other examples: the 3:2 and 5:4 resonances. For the 3:2 resonance, we have 
hve terms dehning (see Table 3), among which Tsao-i, T 3311 , % 3 io are dominant 

in different regions of the (e,z)-plane. In Figure 4 top left, the term 7330-1 of RUrfh 
dominates and using (3.4) we conhrm an amplitude of 7.45 km for e = 0.1, i = 10°. 
For e = 0.1, i = 70° we hnd from Figure 1 that 73311 dominates, while (3.4) yields an 
amplitude of 8.71 km in full agreement with Figure 4, top right. 

Next we analyze the behavior of the 5:4 resonance, where we have hve terms dehning 
Rlarth^ (see Table 3). In the bottom row of Figure 4 we provide the FLI for the 5:4 
resonance as a function of semimajor axis and inclination. Provided that we select regular 
regions, the amplitude of the resonant islands is in good agreement with the size given 
by (3.4). For example, let us hx i = 60° and the eccentricities e = 0.005 and e = 0.5. 
Then, from Figure 1 we infer that the dominant terms are, respectively, Tesio and Tssn. 
Their amplitudes, as computed through (3.4), are about 0.53 and 2.98 km in agreement 
with Figure 4, thus yielding a further conhrmation of the validity of the estimate (3.4), 
when dealing with regular motions exhibiting a pendulum-like structure. 
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4. Detecting the splitting or superposition of resonances 

As mentioned in Section 2, the quantities Snmpq in (2.2) depend on the angle "^nmpq in 
(2.3). The variation of '^nmpq depends on the frequencies cj, D, which can be small, but 
not exactly zero, due to the effect of the secular part^. As a consequence, for a specihc 
resonance, the angles "^nmpq for different n, m, p, q are stationary at different locations. 
As already remarked in [3], this means that each resonance splits into a multiplet of 
resonances. As a consequence, each harmonic term of a specihc resonance, with big 
enough magnitude, yields equilibria located at different distances from the center of the 
Earth. When the width of the resonance associated to each component of the multiplet is 
smaller than the distance separating these resonances then a splitting phenomenon takes 
place, otherwise we have an opposite phenomenon, called superposition, which gives rise 
to very a complex dynamics. 

We also remark that the values provided in Table 2 give just a hint on the location 
of the minor resonances. Indeed, the position of the minor resonances, as well as the 
regular and chaotic behavior of the corresponding resonant regions, are strongly affected 
by the interaction between the secular and resonant parts. A thorough investigation of 
splitting and superposition of resonances is provided in the next section. 

As an example of splitting and superposition of resonances, we consider the 5:3 res¬ 
onance for two different sets of values of the eccentricity and inclination. Besides the 
islands due to T 5510 and 7550-2 in Figure 5, upper left, located at a = 29996.3 km and 
a = 29998.1 km, there appear two thin structures at a = 29997.1 km and a = 29995.5 
km, associated to Tesi-i and 76521, respectively. For larger eccentricities and inclinations 
(Figure 5, upper right) the islands due to 765i-i and 76521 overlap with the main island 
associated to 75510 . 

4.1. An algorithm for distinguishing between splitting and superposition. We 

analyze a specihc resonance for which the dominant terms have been identihed in Sec¬ 
tion 2. For each component of the multiplet we can estimate the corresponding amplitude 
by means of (3.4). Now, we proceed to determine carefully the location of the center of 
the islands, so that the knowledge of the centers and the amplitudes will easily allow 


^Since the coefficient J 2 is much larger than any other zonal harmonic coefficient (see Table 1), the 
secular part is dominated essentially by the J 2 harmonic terms. Without loss of generality, it is enough 
to discuss here just the influence of the J 2 harmonic terms in order to catch the main effects of the 
secular part. 
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US to decide whether we are in presence of a splitting or rather a superposition of the 
resonances. 

For a resonance m : {n — 2p + q), let us write the dominant term Tnmpq in the form 

Tnmpq Cs(^(Tm,n—2p+q 

for a suitable function A = A{L,G, H) and where 

(^m,n-2p+q = {u - 2p + q)M - m9 + {ji - 2p + q)uj + mVt ; 

as in Section 3.1 cs can be either sine or cosine. We look for equilibria satisfying the 
equations: 

L = 0 

^m,n—2p+q 0 . 

Let US consider just the contributions of the secular part and of the dominant term Tnmpq, 
so that we can write the corresponding Hamiltonian i'^ form 

V}-Zl'^'^\L,G,H,am,n- 2 p+q,Uj) = -^+Rsec{L, G, H)+A{L, G, H) cs{am,n- 2 p+q-quj-m\nm) . 
Then, we have that L = 0 if 


^m,n—2p+q b ) 

where 7 = 0 (mod. tt) if cs is cosine and 7 = 7 r /2 (mod. vr) if cs is sine. Equation (4.1) 
determines the equilibria and, in particular, the center of the resonant island. At the 
equilibria we hnd: 


a. 


m,n—2p-\-q 


= {n — 2p + q)M — m + {n — 2p + q)uj + mQ 


= (n — 2p + q)^ — m + (n — 2p + q)( 
± [{n - 2p + q){—+—) + m— 


dL 


+ 


OR.. 


dG 


’ \ dRsec 


dL dG' dH' ’ 

where the ± depends on which equilibrium point we are considering and whether cs is 
sine or cosine. From the condition &m,n- 2 p+q = 0 we compute the value of the semimajor 
axis, which corresponds to the center of the island. 

At this point we have all the ingredients to investigate whether the islands associated 
to two different terms, say R = Tmpq and ^2 = Tn'm'p'q', are splitting or overlapping 
(compare with [ 6 ]). Assuming that the centers of the two islands have coordinates (oi, Ui), 
( 02 , 0 - 2 ) with (Ti = ( 72 , let Ai(e, i), A 2 (e, i) be the amplitudes of the corresponding islands. 
Let D = |ai — 02 ! be the distance between the centers. Then, if (Ai + A 2)/2 < D, we have 
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that the two islands are well separated, while if (Ai + A2)/2 > D the two islands overlap. 
This simple compntation allows us to predict the behavior of neighboring islands. 

4.2. The 4:1 resonance. As an example of the application of this criterion, we consider 
the 4:1 resonance and we £x e = 0.1, ca = 0°, 12 = 0°, while we consider two values of the 
inclination, i = 35° and i = 50°. 

For the 4:1 resonance, the seven terms dehning Rl^rth^ are listed in Table 3. For moder¬ 
ate inclinations, say between 20° and 60°, the term 75420 is dominant for all eccentricities. 
Moreover, excluding the inclination i = 78.5°, this term is also dominant for large incli¬ 
nations, provided the eccentricity is small enough. Table 5 provides the values Uc for the 
centers and the amplitudes A of the dominant terms for i = 35° and i = 50°; moreover 
we report also the distance D from the largest term 75420 - 


Dominant term 

<Tc 

A (km) i = 35° 

D 

A (km) i = 50° 

D 

7441-1 

301.40 or 121.40 

1.17 

3.15 

1.41 

1.42 

74421 

301.40 or 121.40 

1.01 

3.15 

1.80 

1.42 

7541-2 

80.43 or 260.43 

0.216 

6.30 

9.47-10-2 

2.85 

75420 

80.43 or 260.43 

2.73 

- 

3.386 

- 

75432 

80.43 or 260.43 

0.204 

6.30 

0.40 

2.85 

7642-1 

79.66 or 259.66 

1.01 

3.15 

0.38 

1.42 

76431 

79.66 or 259.66 

1.08 

3.15 

1.44 

1.42 


Table 5. Resonance 4:1. Values ac for the centers, amplitudes A in 
kilometers of the dominant terms for i = 35° and i = 50°, distances D 
from the largest term 75420 associated to the 4:1 resonances for e = 0.1, 
a; = 0°, = 0°. 


Although there are seven terms, the number of resonant islands is five (see Figure 5, 
middle left panel) because the arguments of 7441-1 and Tu 2 -i on the one hand, as well as 
74421 and 76431 on the other hand, are the same (modulo a constant). Therefore, 744i_i 
and 7642-1 give rise to a single resonant island with the stable point located between 
(741 = 4 A 44 — a; -|- 360° = 301° — u and a^i = 4X^4 — uj + 180° = 260° — uj. The same 
happens for 74421 and 76431- According to the criterion presented before, it is easy to 
check that all multiplets split for i = 35°, while for i = 50° the term 75420 overlaps with 
7441 - 1 , 74421 , 7642 - 1 , 76431 . This result is validated by the FLI plots provided in Figure 5, 
middle panels. 
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The phenomena of splitting and snperposition of resonances occnr for most of the 
minor resonances stndied in this paper. However, for resonances located at increasingly 
large distances from the Earth the splitting phenomenon becomes less evident, snch as 
the cases of the 4 : 3 and 5 : 4 resonances, or even absent as in the case of the 3 : 2 
resonance. 

When the inclination is eqnal to i = 63.4°, a valne called the critical inclination, 
the argnment of perigee becomes constant (see [14, 3]). Since the argnment of any two 
harmonic terms differs by an integer mnltiple of u, then the shift in semimajor axis is 
zero. As a conseqnence, we conclnde that for the critical inclination (and for very close 
valnes) the pattern of the resonance has a pendnlnm-like strnctnre. 

Using the phenomenon of splitting and snperposition of resonances, we can propose 
a mechanism of transfer from one region to a nearby one by increasing the eccentricity 
or the inclination, and by nsing the snperposition of the islands associated to the differ¬ 
ent dominant terms to move the objects with a minimnm effort. This mechanism conld 
be snccessfnlly applied when the dynamics is like that shown in Fignre 5, middle left 
panel, where there is a coexistence of several nearby distinct islands. However, chang¬ 
ing the orbital plane is dehnitely an expensive manenver (see [7]). A cheaper solntion, 
adopted also in some space missions, consists in modifying the argnment of the perigee 
(see [7], [9]). A change of the argnment of the perigee is shown, for example, in Fignre 5, 
lower panels; the conseqnence of snch change is an evident modihcation of the strnctnre 
of the dynamics of the 4:1 resonance between u = 0° (left plot) and u = 150° (right plot). 


4.3. The 5:1 resonance. As a farther example, we consider the 5:1 resonance for which 
we have hve terms dehning RUrth^ (see Table 3). The terms 7551-2 and 7552-1 prevail 
for small inclinations, otherwise 75520 is dominant. Since 75520 is of order of nnity, while 
the other terms dehning ^^e of order 0{e) or Ci(e^), then for small eccentricities 

pendnlnm-like plots are obtained for each inclination. The stable point is located at 
(751 = 5 A 55 = 105°. 

The 5:1 resonance tnrns ont to be interesting for moderate and large eccentricities, 
where both splitting and snperposition phenomena are clearly distingnished. Thns, for 
small enongh inclinations we have the splitting phenomenon (Fignre 6 , left), while for 
large inclinations the resonances overlap (Fignre 6 , right). 
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For a given eccentricity, there is an inclination where resonances are no longer sep¬ 
arated, but they start to overlap. This inclination can be determined analytically by 
comparing the shift in semimajor axis of the location of the equilibria and the ampli¬ 
tudes of the terms dehning R^arth- More precisely, from Figure 6, left, it follows that 
the islands with the largest width are those associated to 75520 (at o = 14412 km for 
e = 0.2, respectively at a = 14407 km for e = 0.5) and 7652-1 (at a = 14417 km for 
e = 0.2, respectively at a = 14414 km for e = 0.5). Denoting by Ai{e,i) and A 2 {e,i) 
the amplitudes of the resonant islands associated to 75520 and 7652-1, respectively, and 
by D{e,i) the distance (in semimajor axis) between the equilibrium points associated 
to these islands, then, as described in Section 4.1, the superposition takes place when 
D{e, i) < (Ai(e, i) + A 2 (e, i))/2. 

5. Transcritical bifurcations 

The occurrence of transcritical bifurcations is a well known phenomenon which indi¬ 
cates that the stability is transferred from one equilibrium point to another. Transcritical 
bifurcations are very common in almost all minor resonances. 

By analyzing each single term associated to the different minor resonances, one can 
have several examples where the following happens for a given inclination for i < io 
there exist two equilibrium points, one stable and the other unstable; a.t i = io the two 
equilibria annihilate each other; for i > io the stable point becomes unstable, while the 
unstable equilibrium becomes stable. 

Although the theory of transcritical bifurcations is well known, let us make an explicit 
example to clarify how this notion can be applied to minor resonances. Let us simplify 
the discussion by retaining only one term at time in the series development around a 
given resonance. In particular, we consider a Hamiltonian function of the form 

77(S, a) = h{E) + J /(S, i) sin(a - a) , (5.1) 

where d is a constant, 7 is a small parameter (precisely, it coincides with any of the Jnm), 
S is the action conjugated to the resonant angle a, h is the purely gravitational Keplerian 
part, and /(S, i) is a function depending also on the inclination (equivalently, one may 
assume that / is a function of the eccentricity in order to get transcritical bifurcations 
as the eccentricity varies). Then, Hamilton’s equations associated to (5.1) are given by 

S = — 7/(S, i) cos((T — d) 

^ 5/(S,i) . ^ 

= w[^) + 7 ——— sm(cr - a) , 


cr 


(5.2) 
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where w{E) = The stationary points are given by the pairs where ag is 

snch that cos((T — a) = 0, i.e. = a + | or = a + |7r, while is snch that the 
right hand side of the second equation in (5.2) is zero at a = j = 1, 2. 


To look for the linear stability, we compute the eigenvalues 6 of the matrix 



Since cos(cr — d) = 0 at equilibrium, neglecting 0{J‘^) we obtain the following secular 
equation in the variable <5: 



Then, for = d + | we obtain 



while for af'* = d + |7r we obtain 



This shows that if at i = zq fhe function / = /(S,i) reverts sign, then and 


change their stability. The discussion can easily be adapted to the case where the sine 
in (5.1) is replaced by a cosine. 

For each minor resonance, we report in Table 6 the harmonic terms which change their 
sign, together with the inclinations at which this event happens. This does not mean 
that any inclination Iq quoted in Table 6 is automatically a transcritical bifurcation 
point, because we do not know in advance if the harmonic term in question gives rise to 
equilibrium points for inclinations close to 

We expect a bifurcation phenomenon to happen when either the harmonic term that 
changes its sign for i = io is also dominant in some regions located close to Iq and 
moreover all other resonant harmonic terms are small in magnitude in that regions, or 
either the inclination io is such that the splitting phenomenon takes place and the given 
harmonic term is sufficiently large to generate resonant islands for inclinations close to 
if). Since the splitting phenomenon occurs for small inclinations (see Section 4), while the 
inclinations reported in Table 6 are large, the latter case is impossible for all considered 
minor resonances. However, the conditions specified for the former case are satisfied for 
many resonances. Indeed, in Table 6 we report in bold the harmonic terms which are 
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also dominant in some regions of the (e,i)-plane. From these terms, just the underlined 
ones are dominant in some regions close to io (see Figure 1). A detailed analysis shows 
that, within these regions, the underlined terms have a magnitude much larger than any 
other resonant term and give rise to a transcritical bifurcation. As specihc examples we 
present in detail the 4:3 and 5:2 resonances. 



Term 

*0 

Term 

*0 

Term 

*0 

Term 

*0 

3:1 

7431-1 

60° 

74321 

0 

0 

- 

- 

- 

- 

3:2 

74310 

0 

0 

0 

74322 

90° 

- 

- 

- 

- 

4:1 

7541-2 

53.r 

75420 

78.5“ 

7642-1 

51.9'’ and 87.2'’ 

76431 

72.5° 

4:3 

75410 

53.1'’ 

75422 

78.5° 

- 

- 

- 

- 

5:1 

7652-1 

70.5'’ 

76531 

90° 

- 

- 

- 

- 

5:2 

7651-2 

48.2'’ 

76520 

70.5'’ 

76532 

0 

0 

- 

- 

5:3 

7651-1 

48.2° 

76521 

70.5° 

- 

- 

- 

- 

5:4 

76510 

48.2'’ 

76522 

70.5° 

- 

- 

- 

- 


Table 6 . Harmonic terms changing their sign for i = io] those in bold 
are also dominant for some parameter values (see Figure 1). The 
underlined terms are dominant in regions close to io- 


5.1. The 4:3 and 5:2 resonances. As an example, we consider the 4:3 resonance, 
which has hve terms dehning RUlth (^^e Table 3). 

Excluding the inclinations i ~ 0° and i ~ 53.1°, 75410 is dominant for small eccentrici¬ 
ties. For moderate and large eccentricities, we have a balance between two terms, namely 
7440-1 and 74411 (see Figure 1). At i = 53.1° a transcritical bifurcation takes place for 
small eccentricities, as it is shown in Figure 7, top panels. 

For the 5:2 resonance, we have hve terms dehning R^arth (®®® Table 3). The term 
75521 is dominant for large inclinations, provided the eccentricity is large enough, while 
76520 is dominant in the rest of the (e, i)-plane, excluding some small inclinations. For 
small eccentricities, a bifurcation phenomenon takes place at i = 70.5°, as it is shown in 
Figure 7, bottom plots, which provides the FLI values for i = 60° and i = 80°. 

6 . A MORE COMPLETE MODEL 

The purpose of this section is to complement the study realized by using the Hamilton¬ 
ian formulation with some results obtained in Cartesian coordinates on a more complete 
model, which is not limited to the geopotential. 
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We perform a numerical integration in Cartesian variables, including, besides the 
geopotential, the gravitational attraction of Sun, Moon and solar radiation pressure. 
In this way we validate the Hamiltonian model and verify the results obtained in the 
previous sections. 

Concerning the Hamiltonian formulation, we removed in Section 2 the short periodic 
perturbations by averaging over the fast angles. The averaged Hamiltonian contains 
secular and resonant terms, leading to the determination of the mean orbital elements. 
Therefore, for the equations of motion in Cartesian coordinates, in order to represent the 
FLI as a function of the same variables, we transform from osculating orbital elements to 
mean elements. This computation implies a numerical average of the osculating elements, 
which is performed in the course of the integration itself. 

We stress that each of the disturbing forces due to the geopotential. Moon, Sun and 
solar radiation pressure induces a short periodic variation of the orbital elements. The 
stronger effects are notably due to J 2 , since the short periodic harmonic terms of order 
J 2 are much larger in magnitude than any other short periodic term. 

The results obtained by using the Hamiltonian formulation are validated by integrating 
the Cartesian equations of motion as in Figure 8 . We remark that the computation of 
Figure 8 takes a machine execution time 12 times longer than the plot obtained using the 
Hamiltonian formalism. We select the following resonances: 3:1, 3:2, 5:1, 5:3. We have 
used as starter a single step method (a Butcher numerical algorithm), while a multistep 
numerical method (Adams-Bashforth 12 steps and Adams-Moulton 11 steps) performs 
most of the propagation. All hgures are obtained for a dynamical model which includes 
also the gravitational attraction of Sun, Moon and solar radiation pressure. For the 3:1 
and 3:2 resonances we considered the Earth’s gravitational potential up to degree and 
order n = m = 4. For the 5:1 and 5:3 resonances, we considered also the effects of J 55 
and Jes. 

The panels of Figures 8 must be compared with those obtained using the Hamiltonian 
approach, precisely Figures 3 (middle row, left panel), 4 (upper left panel), 6 (right 
panel), 5 (upper left panel). 

The comparison leads to the following conclusions: all dynamical features of the minor 
resonances, which were explained by using the Hamiltonian formalism, are retrieved by 
integrating the full equations of motions; the perturbations due to Sun, Moon and solar 
radiation pressure with a small A/m parameter do not modify signihcantly the main 
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characteristics, like the location of the equilibrium points, the amplitude of the resonant 
islands and the regular or chaotic behavior of the orbits. 


Appendix A. Secular and resonant terms 

We report below the explicit expressions of the terms which provide the secular part 
and the resonant parts appearing in Table 3. 

The leading terms of the expansion of the secular part are the following (notice that 
the hrst two terms are zero, since they are of third order in the eccentricity): 


7200-2 — 72022 — 0 , 
I^eR%J2 ( 3 


72010 — 






_ p2'|-3/2 

3 


sin^i —-sinije(l — e^) ^^^sino;. 


Tsoi-l — 73021 — 

^ ^ I-1‘ER%J4: ( 35 . 4 . 15 .2 A /I 2'i-7/2 /r) \ 

7101-2 = 71o32 = — sm 7 + — sm zj — (1 - e ) cos(2n;), 

^ieR^Ja /105 . 4 . 15 . 2 • I 3\, 2\-7i2 

74020 =-;- -sm z-sm z + - 1 H-1 — e ' . 

“ “ a5 V 64 8 8/^2^^ ’ 


Let (Tsi = M — 30 + a; + 311; we report below the leading terms of the expansion of the 
3:1 resonance: 


7330-2 = 
73310 = 
73322 = 
7431-1 = 
74321 = 


^j^eReJss { 15 


fiER^JsS f 

45 . , 

1 

8 

^IeRe^SS f 

45 . , 

1 

8 

IJ^eR'eJas 

f 105 

■ 

1 8 " 

I^eRe'^^s 

1-315 


1 8 


|y(l + cosz)^ — COs((T3i +2u- 3A33)| , 


lle^ 
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Let (T 32 = 2M — 36* + 2a; + 3fi, we report below the leading terms of the expansion of 
the 3:2 resonance: 


Tsso-i — 


r 15 


— (1 + cosi)'^ ecos((T 32 + a; - 3 A 33 ) \, 


Tsaii — 


f 45 . 


sin^i(l + cosi) 3 ecos(cr 32 — a; — 3A33)|, 

J43 /105 . , ..3 ox 

/430-2 =-?-f — sin *(1 + cosz) — sin((T32 + 2a; - 3A43) >, 

I 16 2 J 

sini(l — 3 cos^i — 2cos^i)(l + e^) sin((T32 — 3A43)|, 


74310 — 


a5 I 

?4 


^eReJ^'a r 315 .3. 2 ■ ( o Q^^\ 

/4322 =-R-i —^ Sin'" I cost be sin((T32 - 2a; - 3A43) > . 

I 8 ) 


Let (741 = M — 40 + a; + 411; we report below the leading terms of the expansion of the 
4:1 resonance: 


7441-1 — 
74421 = 


P.E-Re744 I 1^ ^.^2 (. Qg ^^2 6 cos(cr4i + 02 — 4A44)| , 


5 i 4 

I^^eR-e^^^ f 315 . ^ , 5e 


O'" 

?4 


I — sin"^ i y cos((T 4 i - o; - 4A44) | , 


1^eReJ5a (2835 


7541-2 = ——a - 1 -sin i (3 + 4 cos i — 6 cos^ z — 12 cos'^ z — 5 cos"^ z) 

I 256 

X sin((T4i + 2a; — 4 A 54 ) [, 


75420 — 


75432 — 


/Ze 72|;754 r 945 , . 


ry4o „ 0 . . Q, 4-\ /-I I 06 \ 

< —— sm z (1 — 4 cos z — 6 cos z + 4 cos z + 5 cos z) 1 H-— 1 

I 16 V 2 / 

X sin( 0-41 - 4 A 54 )| , 

/^E-R|;754 J" 27405 •,« 2-|3 3- e; 4 -12 

-j-<-smz ( — 1 — 4cosz + 6 cos z + 4cos z — 5cos z) e 

a 6 I 128 ^ ^ 

X sin((T 4 i — 2a; — 4 A 54 ) [, 


7642-1 — 


^ieR^Jga ( 3e 


945 


(—1 — 2 cosz + 2 cos^ z + cos^ z) + 


10395 

128 


(1 + 12 cosz 


+ 6 cos"^ z — 20 cos'^ z — 15 cos"^ z) sin^ z cos(cr 4 i + a; — 4A64) \ , 


76431 — 


^ieR%Jm ( 7e 


945 


10395 


(—3 + 6 cos^ z — 3 cos"^z) + 

2 L 16 ^ ^32 


(1 — 6 cos^ z 


+ 5 cos'^ z) sin^ z cos(cr 4 i — cu — 4A64) 
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Let (743 = 3M — 46* + 3a; + 4fi; we report below the leading terms of the expansion of 
the 4:3 resonance: 

4 

75 


^eReJaa f 105 •^4 3e f 4\ 

7440-1 =-y C0 s((743 + O; - 4 A 44 ) J , 




^eR%Jaa r 105 . 2 I •\2 Oe / /i\ 

74411 = --I ~ ^ y COS(cr43 -OJ- 4A44) I , 


7540-2 — 


^eR%J^A [8505 . . ^ •^4 2 • ^ I o 4 1 ll 

-it;- \ * (1 + *^ 0 ^ *) ^ Sin(cr 43 + 2 a; - 4 A 54 ) \, 


^lER%JbA [945 ■■ lo i A ■ r 2 - 10 3- k 4•^ 

/5410 =-gf Sin ^ (3 + 4 cos z — 0 cos z — 1/ cos z — 5 cos z) 

CL C oZi 

X (^1 - — j sin(a 43 - 4 A 54 ) | , 

llER%JbA [82215 ■ ■ ^ ■ p 2 - 1/1 3 -ic; 4•^2 

75422 =-^f- Sin z 11 — 4 cos z — 0 cos z + 4 cos z + 5 cos z I e 

a® I 128 ^ ’ 

X sin ((743 — 2a; — 4 A 54 ) 

Let (751 = M — 50 + a; + 5f2; we report below the leading terms of the expansion of the 
5:1 resonance: 

fiER^Jbb ( 14175 2 /-, Q -lO 2-0 3-0 4 - 5 -\ 

7551-2 =-g- 5 -e (1 + 3 cos z + 2 cos z — 2 cos z — 3 cos z — cos z j 

CL t 256 


75520 — 


X cos((75i + 2a; - 5A55)|, 

f^ER%J55 [9450^ 13e^/^^ -o 2-0 3-1 a ■ , 5 ■\ 

-^-<- 1 H-1 + cosz — 2cos z — 2cos z + cos z + cos z 

a6 I 32 ^ 2 ^ 

X cos(c75i - 5 A 55 ) [ , 


75532 — 


/z_b-R|; J 55 [ 274050 2 


1- 


e^(l — cos z — 2 cos^ z + 2 cos^ z + cos'^ z — cos^ z) 


256 


X cos(( 75 i — 2a; — 5 A 55 ) [, 


7652-1 — 


^^ER%J(i5 [ 155925 


e (1 — cos z — 6 cos^ z — 2 cos'^ z + 5 cos"^ z + 3 cos“ z) 


76531 — 


I 128 

X sin z sin(( 75 i + a; — 5A65) \ , 
fJ^ER%J65 [ 1455300 


V 


-e(— cos z + 2 cos'* z — cos^ z) 
128 ^ ^ 

X sin z sin((75i — u — 5 A 65 ) 


Let (752 = 2M — 50 + 2a; + 5f2; we report below the leading terms of the expansion 
of the 5:2 resonance (notice that the hrst term is zero, since it is of third order in the 
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eccentricity): 


Tssi-i — 0, 

^lER%Jbb (9450 


75521 — 


^-<-e(l + cosi — 2cos i — 2cos'^i + cosS + cos^i) 

6 I 8 

X cos(a52 -u- 5 A 55 ) \ , 


( 10395 2 / 
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^eR%Jg5 J 51975 13e^ . 2 ■ ^ 3-,r 4 ■. o 5-', 

76520 =-^ i :—(1 H-)(1 — cosz — Dcos * —2cos * + 5cos z + 3cos t) 

a' I 64 2 

X sinisin((T52 - 5A65)| , 


76532 — 


fJ^ERE^Gb f 3638250 2/- • I o 3 • 5 

- 7 - 1 — 75 E—e (—cos* + 2cos * — cos *) 

ci‘ I 128 

X sin*sin((T 52 — 2a; — 5 A 65 ) [ . 


Let (T 53 = 3M — 50 + 3a; + 5f2; we report below the leading terms of the expansion of 
the 5:3 resonance: 


7550-2 — 


j^ER%Jbb 18505 2^1 I •^5 i I o e;\ 

- 6 -1 ”577“^ (1 + COS *) cos((T 53 + 2a; - 5 A 55 ) \, 




fJ^ER^Jbb f4725^ IQ -lO 2- 0 3- Q 4- 5•^ 

75510 = - Q -< (1-—)(1 + 3cos* + 2cos * —2cos * —3cos * — cos *) 

CL V ^ 


75522 — 


X COs((T53 - 5A55)| , 

f^ER%Jbb [822150 27t I • o 2 • 0 3 • I 4 • I 5-1 

a6 I 256 ^ ^ 


X cos((T 53 - 2a; - 5 A 55 ) \, 


7651-1 — 


/4_e-R|;J 65 r 10395 


<--—e (2 + 5 cos * — 10 cos^ * — 10 cos"^ * — 3 cos^ *) 

I 64 


76521 — 


X sin * sin((T 53 + a; — 5A65) |, 

/^i?-7?^765 [571725 , . 2 ■ r) s-ik 4-iq 5-'i 

' 6(1 — COS * — 6 cos * — 2 cos * + 5 cos * + 3 cos *) 


1 - 


128 


X sin * sin((T 53 — cu — 5 A 65 ) 
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Let (754 = 4M — 56* + 4a; + 5fi; we report below the leading terms of the expansion of 

the 5:4 resonance: 

/iE-R|J 55 / 94 5 5 e:\ 

7550-1 =--1 ^ cos(a 54 + o; - 5 A 55 ) , 

j 14175 /-. Q -lO 2 - o 3 - Q 4 - 5 -', 

/ 55 II = - - -< -e( 1 + 6 cos I + Z cos I — Z cos i — 6 cos I — cos l) 

a° I 16 

X COs((754 -Ul - 5A55)| , 

llER%JGb r 10395 2 /I I -^5 • • • ^ I o e;\ 

7650-2 =--1 e (1 + cosz) sin*sin(a 54 + 2a; - 5A65)| , 

Tesio =- j —Qo“ = —o”^ ) (2 + 5 cos ? — 10 cos ? — 10 cos i — c> cos i) 

CL V ^ 

X sinisin((754 — 5A65)|, 

fiER^Jbb 1987525 2^-1 . ^ 2 • r) 3 • I K 4 • I Q 5 

/6522 = - ^i -e (1 — cos t — b cos t — Z cos z + 5 cos z + 3 cos z) 

o' I 64 

X sinzsin ((754 — 2a; — 5A65)| . 
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Figure 3. FLI for the 3:1 resonance for w = 0°, = 0°: i = 10° in the 

top row and i = 30° in the middle row, e = 0.005 in the left column and 
e = 0.5 in the right column. In the bottom row we provide the FLI as a 
function of inclination and semimajor axis for <731 = 243°, cu = 0°, 12 = 0° 
for e = 0.005 (left) and e = 0.5 (right). 




















28 


A. CELLETTI AND C. GALE§ 


32195 

32190 

32185 

32180 

32175 

32170 

32165 

32160 



32195 

32190 

32185 

32180 

32175 

32170 

32165 

32160 



50 100 150 200 250 300 350 400 


8.5 
8 

7.5 
7 

6.5 
6 

5.5 
5 

4.5 
4 

3.5 



Figure 4. Upper row: FLI for the 3:2 resonance for a; = 0'’, = 0'’, 

i = 10°, e = 0.1 (left); i = 70°, e = 0.1 (right). Bottom row: FLI for the 
5:4 resonance as a function of inclination and semimajor axis for uj = 0°, 
Vt = 0°, (754 = 105°, e = 0.005 (right); e = 0.5 (right). 
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Figure 5. Upper plots: FLI for the 5:3 resonance with a; = 0°, hi = 0°, 
e = 0.1, i = 15° (left); e = 0.3, i = 45° (right). Middle plots: splitting and 
superposition of resonances. FLI for the 4:1 resonance for e = 0.1, cj = 0°, 
n = 0°: i = 35° (left); i = 50° (right). Lower plots: FLI for the 4:1 
resonance for e = 0.1, i = 42.5°, hi = 0°: a; = 0° (left); u = 150° (right). 
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Figure 6 . Upper plots: FLI for the 5:1 resonance for cj = 0'’, hi = 0'’: 
i = 30°, e = 0.2 (left); i = 45°, e = 0.2 (right). 
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Figure 7. Upper panels: FLI for the 4:3 resonance for e = 0.005, 
cn = 0°, n = 0°: i = 40° (left); i = 70° (right). Lower panels: FLI for the 
5:2 resonance for e = 0.005, cn = 0°, hi = 0°: i = 60° (left); i = 80° (right). 
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Figure 8. FLI (using Cartesian equations) under the effects of the 
geopotential + Moon + Sun + SRP with A/m = 0.01 [m^/%]. Top left: 
the 3:1 resonance for e = 0.005, i = 30°, cn = 0°, hi = 0°; Top right: the 
3:2 resonance for e = 0.1, i = 10°, a; = 0°, hi = 0°; Bottom left: the 5:1 
resonance for e = 0.2, i = 45°, a; = 0°, hi = 0°; Bottom right: the 5:3 
resonance for e = 0.1, i = 15°, a; = 0°, hi = 0°. 












